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Abstract. An optimization of the CERN SPL beam line has been performed guided by the sensitivities 
to the #13 mixing angle and to the Sep Dirac CP violating phase. A UNO-like 440 ktons water Cerenkov 
detector located at 130 km from the target in a new foreseen Frejus laboratory has been used as a generic 
detector. Concerning the Sep independent 0i3 sensitivity, a gain of about 20% may be reached using a 
3.5 GeV proton beam with a 40 m long, 2 m radius decay tunnel compared to the up to now consid- 
ered 2.2 GeV beam energy and 20 m long, 1 m radius decay tunnel. This may motivate new machine 
developments to upgrade the nominal SPL proton beam energy. 



1 Introduction 

The very near future of the neutrino long baseline experi- 
ments is devoted to the study of the oscillation mechanism 
in the range of Am 2 = Am 2 atm w 2.4 x 10~ 3 eV 2 [HE] us- 
ing conventional beams. The current K2K experiment 
in Japan [2], and the forthcoming MINOS in the USA 
jH] take benefit of low energy beam to measure the Am 2 
parameter using the disappearance mode — > i/^, while 
OPERA/ICARUS experiments |4ll5] using the high energy 
CNGS beam [g] will be able to detect v T appearance. If 
we do not consider the LSND anomaly [2 that will be 
further studied soon by MiniBooNE experiment (Sj, the 
three flavor family scenario will be confirmed and accom- 
modated by a 3 x 3 Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix jg] with three angles (6*12,^13,^23) 
and one Dirac CP phase Sep- 

Beyond this medium term plan, two of the next fu- 
ture tasks of neutrino physics are to improve the sensi- 
tivity of the last unknown mixing angle parameter, the 
so-called #13, and to explore the CP violation mechanism 
in the leptonic sector. The present upper bound on 613 is 
sin 2 26»i3 < 0.14 for Am 2 = Am 2 atm (90% CL) pE]. This 
sensitivity can be improved using reactor and accelera- 
tor experiments. In reactor experiments, one uses T> e in 
disappearance mode and may reach sin 2 26*13 < 0.024 for 
Am 2 — Am 2 tm (90%CL) [U]. In accelerator experiments, 
and D e from (3 beams [12] in both disap- 
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beams either with sub-mega watt proton drivers [131114] or 
with multi-mega watt proton drivers |14U15Ulfi| . The later 
neutrino beam type, called Superbeam, is foreseen to be 
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extended to produce beam and beam from muon 
decays, the so-called Neutrino Factory, in order to study 
the eventual leptonic CP violation. Such neutrino complex 
is under study in Japan, in USA and also in Europe at 
CERN and details may be found in reference ^H] . A com- 
parison of the performances of (3 beam and Superbeam 
may be found for instance in reference 23 • The reactor 
experiment result on 6*13 is straight forward as compared 
to Superbeam and Neutrino Factory results that are on 
one hand richer but in an other hand more complex to 
analyse due to the interplay between the different physics 
factors 6»i3, S C p, sign(/Amf 3 ), sign(tan(202 3 )) [T7HT%j . 

This paper presents results of a new simulation of the 
SPL (Super Proton Linac) Superbeam that could take 
place at CERN [V^, using for definitiveness a UNO-like 
440kT fiducial water Cerenkov detector 20 located in a 
new enlarged underground laboratory under study in the 
Frejus tunnel, 130 km away from CERN (23]. The SPL 
neutrino beam is created by decays of pions, muons and 
kaons produced by the interactions of a 4 MW proton 
beam impinging a liquid mercury jet JH]. Pions, muons 
and kaons are collected using two concentric electromag- 
netic lenses (horns), the inner one and the outer one are 
hereafter called "Horn" and "Reflector" respectively (22]. 
The horns are followed by a decay tunnel where most of 
the neutrinos are produced. A sketch of the beam line is 
shown on figure Q 

The analysis chain consists of different stages: the sim- 
ulation of the interactions between the proton beam and 
the mercury target, the propagation of the resulting sec- 
ondary particles through the magnetic field and the mate- 
rials of the horns, the tracking of ^ , K ±fi and until 
they decay, the computation of the neutrino flux at the 
detector site, and finally the statistical analysis. A part of 



2 



Jean Eric Campagne, Antoine Cazes: The #13 and Sep sensitivities of the SPL-Frejus project revisited 



H" linac, 2.2GeV, 4MW Accumulating Ring 




Fig. 1. Sketch of the SPL neutrino Superbeam from CERN to 
the Frejus tunnel. 



Table 1. Liquid mercury jet parameters. 
Hg target 

Hg jet speed 20 m/s 
density 13.546 
Length, radius 30 cm, 7.5 mm 



the simulation chain has already been described in refer- 
ence |28ll24| . 

Compared to recent papers on the same subject |171 
25 . 26], we have reoptimized the Horn and Reflector shapes 
|27j , and introduced the kaon background simulation which 
allows us to update the SPL beam energy. The organi- 
zation of this document follows the simulation chain: the 
interaction between the proton beam and the mercury tar- 
get is presented in the second section. The kaon produc- 
tion is detailed in the third section. The simulation of the 
horns is described in the fourth section, while the algo- 
rithms used to compute the neutrino fluxes are explained 
in the fifth section. Then, the sensitivities to #13 and Sep 
are revisited with new studies about the optimization of 
the proton beam energy, the pion collection and the decay 
tunnel geometry. 



2 Target simulation 

Since hadronic processes are crucial to describe the in- 
teractions of the proton beam on the target, the FLUKA 
simulator [2Sj has been chosen for this first step of the 
simulation. The target used in the present study is a mer- 
cury liquid jet QS] simulated by a cylinder 30 cm long 
(representing two hadronic lengths) and 1.5 cm diameter 
(see table nj. Other types of target are under study jlrij . 
The pencil like simulated proton beam is composed of 10 6 
mono-energetic protons. The beam axis is also the symme- 
try axis of the target and the horns and the decay tunnel. 
Simulations have been performed for 2.2 GeV proton ki- 
netic energy, the up to now nominal design fTJ)], as well 
as for 3.5 GeV, 4.5 GeV, 6.5 GeV and 8 GeV according to 
possible new designs |29] . 



Particle production yields are summarized in table El 
The pion momentum spectra obtained at different energies 
and normalized to a 4MW SPL beam power are presented 
in figure Eta). At low energy, pions come from A decays 
while the high energy part is dominated with multi pion 
production. At very low energy, for P < 200 MeV/c, pions 
come from A produced by protons of the target excited 
by the beam interactions, while for higher energy, pion 
production is due to transformation of protons of the beam 
into A. 

The horns are designed to focus the 600 MeV/c pions 
(see section 0} and the variation of the number of such 
pion is rather smooth with respect to the beam energy 
considering a 4MW fixed beam power: 4.19 x 10 13 7r/s for 
the 2.2 GeV beam, 4.91 x 10 13 tt/s for the 3.5 GeV beam, 
5.14 x 10 13 tt/s for the 4.5 GeV beam, and 4.92 x 10 13 tt/s 
for the 6.5 GeV beam. The main difference is made by the 
angular distribution. Figure Gib) shows this distribution 
for the 7r + exiting the target with a momentum between 
500 MeV/c and 700 MeV/c. The acceptance of the horns is 
limited to the pion below 25°, and we see that more pions 
are accepted by the horns for the 3.5 GeV and 4.5 GeV 
proton beams compared to other beam energies. 

The secondary proton and neutron rates induce impor- 
tant radiation damages and power dissipation in the horns 
which have been addressed in reference |2B], and which 
will require specific R&D effort. At 2.2 GeV, kaon yields 
are very low, but it has a dramatic energy dependence as 
further studied in section E| It is worth mentioning that 
the numbers in table El are not to be taken as face values, 
because the cross sections of pion and kaon productions 
using proton beam are still under studies as for instance by 
the HARP experiment [HOj- The cross section uncertain- 
ties are the main source of discrepancy between simulator 
programs. Some comparisons between FLUKA and MARS 
[32] have already been presented in the same context [23J. 
The energy distribution of the pions exiting the target, 
computed with the two simulator programs FLUKA and 
MARS, is shown on figure Eta). The discrepancy is quite 
large for the low energy part. However, the horns are de- 
signed to focus the high energy part of the spectrum (see 
section QJ, and therefore, MARS and FLUKA are in bet- 
ter agreement for the energy spectrum computed at the 
entrance of the decay tunnel, as shows figureEtb). So, the 
discrepancy at low energy between MARS and FLUKA 
does not matter too much for the present application. A 
difference of 10% has been found between the #13 sensi- 
tivity computed with the two generators (see section 0, 
that can be taken as systematic error. 



3 Kaon production 

The possibility to increase the SPL energy in order to 
study the optimization of the physics program has been re- 
cently pointed out [22j. Then, the kaon production should 
be clearly addressed because it is a source of v e and v e 
background events. The kaon decay c hanne ls and branch- 
ing ratios are presented in table ITU in IA.3I 
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P (GeV/c) 8 (degree) 

Fig. 2. (a) 7r + momentum distribution per second at the exit of the target for the different proton beam energies studied, 
simulated with FLUKA, and (b) ir + angle with respect to the beam axis of the pion having a momentum between 0.5 GeV/c 

and 0.7 GeV/c. The different SPL beam kinetic energies presented are ( ) 2.2 GeV, ( ) 3.5 GeV, ( ) 4.5 GeV and 

( ) 6.5 GeV. 




Fig. 3. 7r + momentum distribution at the exit of the target (a) and at the exit of the horns (b), simulated by FLUKA ( ) 

and by MARS ( ). 
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Table 2. Average numbers of the most relevant secondary particles exiting the 30 cm long, 1.5 cm diameter mercury target 
per incident proton (FLUKA). The fi + / /x~ numbers and the K + /K° numbers have been multiplied by 10 4 . Note that the K~ 
production rate is at the level of 10 -5 per incident proton. 



E k (GeV) 


P 


n 


7 


e+ 


e 


7V+ 


7T 






K+ 


K° 


2.2 


1.4 


17 


5.0 


0.08 


0.17 


0.24 


0.18 


4 


1 


7 


6 


3.5 


1.8 


23 


7.0 


0.15 


0.28 


0.41 


0.37 


10 


3 


35 


30 


4.5 


2.3 


25 


7.7 


0.21 


0.35 


0.57 


0.39 


11 


3.3 


93 


68 


8 


3.1 


33 


11.0 


0.41 


0.63 


1.00 


0.85 


30 


9.5 


413 


340 




Fig. 4. Kaon production (a) as a function of the incident proton beam kinetic energy (Ek) for 10 6 incident protons with ( ) 

curve for K + , ( ) curve for K~~ and ( ) curve for K° . Pion production (b) in the same conditions with ( ) curve 

for 7r + and ( ) curve for ir~ . 



The target simulation described in section has been 
used with 10 6 p.o.t with kinetic energy uniformly dis- 
tributed between 2.2 GeV and 5 GeV. The momenta of 
outgoing pions and kaons are recorded when they exit the 
target. The number of produced K"^ at different proton 
beam energies are presented on figure a). On the one 
hand the K° production rate is similar to the K + produc- 
tion rate, but on the other hand the K~ production rate is 
almost forty times smaller. In comparison, the numbers of 
7r + and ir~~ produced in the same conditions are presented 
on figureBJb). Pion production rate is about two orders of 
magnitude greater than the kaon production rate. The be- 
havior of the two pion and kaon production rates are quite 
different. The ir + yield grows smoothly with the proton 
energy while the production of kaons seems to have two 
origins, which has been confirmed by FLUKA's authors 
[33 . For beam energy below approximatively 4 GeV, the 
resonance production model is used, and one notices a 
low production rate with a maximum at about 3.4 GeV. 
For beam energy above 4 GeV, the dual parton model is 



used, and the production rate experiences a threshold ef- 
fect with a rapid rise. The ratio between positive kaon and 
pion production rates is about 0.5% between 2.2 GeV and 
4 GeV and grow up to 2.3% at 5 GeV. One notices that 
the transition between the two kaon production models 
may not be optimal. 



4 Horns simulation 

The simulation code of the electromagnetic horns is writ- 
ten using GEANT 3.2.1 [33| for convenience and since elec- 
tromagnetic processes are dominant, FLUKA has not been 
considered as mandatory, but this may be revised in a fu- 
ture work. The geometry of the horns has been inspired 
by an existing CERN prototype and a Reflector design 
proposed in reference Depending on the current in- 
jection, only positive secondary particles or negative sec- 
ondary particles are focused. The relevant parameters are 
detailed in table El 
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Fig. 5. Design of the Horn and the Reflector conductors im- 
plemented in the GEANT simulation in case of the generation 
of a 350 MeV neutrino beam. The Hg target is located inside 
the cylindrical part of the Horn. 

Table 3. Relevant parameters of the horns in case of the gen- 
eration of a 260 MeV neutrino beam (or 350 MeV in paren- 
thesis). The shapes of the conductors are independent of the 
proton beam energy, as the focusing has been optimized for a 
600 MeV/c (or 800 MeV/c) pion momentum. 





Horn 


Reflector 


neck inner radius 


3.7 cm 


20.3 cm 


neck length 


40 cm 


120(140) cm 


end cone inner radius 


16 cm 


35.7 cm 


outer radius 


20.3 cm 


40 cm 


total length 


120(140) cm 


190(220) cm 


Alu thickness 


3 mm 


3 mm 


Peak current 


300 kA 


600 kA 


Frequency 


50 Hz 


50 Hz 



The mercury target is localized inside the Horn be- 
cause of the low energy and the large emittance of the 
secondary pions produced: 

< PnT > I < P-k >~ 240 MeV/400 MeV 

(2.2 GeV proton beam energy). This explains the Horn 
design (figure^, with a cylindrical part around the target, 
called the neck, which is larger than the transversal size 
of the target to simulate the room for target handling, 
and a conic part designed such that the relevant pions 
are focused as much as possible to exit the magnetic field 
parallel to the beam axis. 

The shape of the horn conductors is a crucial point 
since it determines the energy spectrum of the neutrino 
at the detector site. The details of the conductor shape 
optimization for the present context may be found in ref- 
erence [22j. We just recall here some ingredients. For a #13 
driven — > v e oscillation, a Am^ parameter value of 
2.5 x 10 _3 eV , and a baseline distance of 130 km, the first 
oscillation probability maximum occurs for a neutrino en- 
ergy of 260 MeV. The optimization of the physics potential 
depends at first approximation on the pion neutrino char- 
acteristics, which energy is fully determined by the pion 
2-body decay and boost. To reach an energy of 260 MeV, 
the pion needs a (3 = 0.97, which in turn induces a pion 
momentum of 600 MeV/c. Then, the shape of the conic 



part of the horns is determined such that these 600 MeV/c 
pions exit parallel to the beam axis. 

An other shape of the horn conductors has been used 
to produce a 350 MeV neutrino beam to compare the sen- 
sitivity potential (see section^ • In that case, keeping the 
current intensity unchanged (300/600 kA), the lengths of 
the Horn and the Reflector should be increased by 16% 
and 18.5%, respectively (see table |2J). 

Before closing this section, it is worth quoting that the 
Horn/Reflector conductor shapes optimized in the present 
study to focus a given pion momentum value, is not af- 
fected at first order by a proton beam energy change. 
What is affected is the production rate of the relevant 
pions. This Horn/Reflector design consideration would be 
different if one wished to focus as much as possible all the 
pions produced for which the mean energy is of course 
affected by a proton beam energy change. 



5 Particle decay treatment and flux 
calculation 

The decay tunnel representation is a simple cylinder with 
variable length (Lt) and radius {Rt) filled with "vacuum" 
and located right after the horns. The default design is a 
20 m long and 1 m radius cylinder, but simulations have 
also been conducted with lengths of 10 m, 40 m and 60 m, 
and radius of 1.5 m and 2 m in the spirit of reference [36J. 
In the GEANT simulation, to gain in CPU time, only 
pions, muons and kaons are tracked in the volume of the 
tunnel, and all particles exiting this volume are discarded. 

Beyond the 1/L 2 solid angle factor due to the source- 
detector distance (L) which decreases dramatically the 
fluxes, the neutrino beam focusing is very limited due to 
the small pion boost factor (« 4). Therefore, computa- 
tional algorithms have been used to avoid a too prohibitive 
CPU time resulting from the simulation of each secondary 
particle decay. Otherwise, about 10 15 p.o.t would have 
been necessary to obtain reliable statistics for the esti- 
mation of the v e flux for instance. 

It is worth pointing out that the particle decays occur- 
ring before the entrance of the decay tunnel are also taken 
into account and treated in the same manner, which is not 
the case in reference [3fi| . 



5.1 Algorithm description 

The decay code has been included in the GEANT code. 
The basic idea of this algorithm is to compute the neutrino 
fluxes using the probability of reaching the detector for 
each neutrino produced by a it or a K or a fi particle (on- 
axis neutrino beam). This method has already been used 
in reference (22] and has been modified and extended to 
the kaon decay chain for the present study. 

Muon neutrino comes mostly from pion decay. In a 
first stage, each pion is tracked by GEANT until it decays. 
Then, the probability for the produced muon neutrino to 
reach the detector is computed. The flux is obtained ap- 
plying the probability as a weight for each neutrino. All 
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Table 4. Number of protons on target for different beam en- 
ergy at 4 MW constant power. One year is defined as 10 7 s. 



Beam energy 
(GeV) 


Number of proton 
per year (10 23 p.o.t/y) 


2.2 


1.10 


3.5 


0.70 


4.5 


0.56 


6.5 


0.40 


8.0 


0.30 



the pions produced in the simulation are therefore useful 
to compute the flux, and this allows to reduce the number 
of events in the simulation to 10 6 p.o.t. In this computa- 
tion, the decay region (horns and tunnel) is considered as 
point like compared to the source-detector distance. 

The same method is applied for neutrino coming from 
muons and kaons with some modifications because most 
of the muons do not decay, and there are very few kaons 
produced (see table [2j). The probability computation is 
presented in appendix El 



5.2 Validation of the algorithm 

The validity of the method presented in the previous sec- 
tion has been tested against a straight forward algorithm 
consisting in decaying each pion N times (N « 10 6 ) in 
a full GEANT simulation of the event (decays included) . 
Such method presents the advantage to keep all the infor- 
mation of the neutrino available for further studies. It can 
be a good approach to compute the muon neutrino flux 
coming from pion decays. It can also provide the beam 
profile, but it shows its limits for the muon induced fluxes, 
especially the v e flux. Indeed, this means that each muon 
is duplicated N times and when a muon decays, it must 
decay N times again. For N sw 10 6 , this is a prohibitive 
CPU time consuming. The and fluxes are displayed 
on figure El for both methods. The two spectra show a 
clear agreement, and this makes reliable the probability 
method. 



5.3 Simulated fluxes 

The fluxes are computed at a distance of 100 km from the 
source by convention and can be rescaled at any desired 
distance. They provide the number of the four neutrino 
species (f M , 9^, v e , v e ) passing through a 100 m 2 fiducial 
area during 1 year. 

In practice, the fluxes are given as a function of the 
neutrino energy via histograms composed of 20 MeV bin 
width. These histograms are filled with the energy of each 
neutrino weighted by the probability to reach the detec- 
tor (section l5.1|l . To obtain the fluxes, the histograms are 
rescaled to the number of p.o.t per year depending on the 
beam energy. Table Q] reports on the number of p.o.t per 
year for the different energies studied using the definition 



of one year being 10 7 s and keeping the beam power con- 
stant (i.e. 4 MW). 

Three origins are identified in the composition of each 
neutrino flux: 

- neutrinos from pions, which includes neutrinos created 

by primary pion decays and neutrinos coming from 
the muons produced by pion decays or muons directly 
exiting the target. This is the component studied in 
reference jHE] but with different settings and event gen- 
erator; 

- neutrinos emitted during the decay chain of the charged 

kaons, either by direct production, or produced by the 
daughter pions and muons; 

- neutrinos coming from the decay chain of the neutral 

kaons. 

The three components of the fluxes for the four neu- 
trino species are presented on figuredfor positive particle 
focusing and a proton beam kinetic energy of 2.2 GeV. 
The flux is dominated by the neutrinos of pion decays, 
but a tail above 500 MeV (insert on the top left part) is 
created by the K + — > [i^v^ channel, which is anyway at 
least three order of magnitude below the flux maximum. 
The Pfj, flux is mostly due to the decays of n~ that are 
not unfocused by the horns, but the higher energy part 
comes from /i + decays. It is noticeable that the v e and v e 
fluxes are respectively more than 200 and more than 7000 
times smaller than the flux. The v e are produced in a 
large part by the K { [ —* ir + e~v e decay channel and by yT 
decays, while the v e flux is dominated by the /i + decays. 

On figure |H1 the horns are set to focus negative parti- 
cles keeping other parameters identical. Comparing with 
positive focusing, one can at first approximation translate 
the results by exchanging particles and anti-particles, ex- 
cept that the K + /K~ ratio is about 50 in the beam-target 
interactions (see table EJ- 

On figures l§l and II II one observes the evolution of fig- 
ured when the proton beam kinetic energy increases to 
3.5 GeV and 8 GeV, respectively. Correspondingly, the re- 
sults for negative particle focusing are presented on fig- 
uresEilandEl One clearly notices the increase of the kaon 
induced neutrino contents as the beam energy grows. 

On tableware reported the integral of the fluxes when 
one modifies the decay tunnel length and radius, as well as 
the beam kinetic energy. Changing the length from 10 m 
to 40 m will increase the flux by 50% to 70% and in 
the same time, the number of v e will be multiplied by 
a factor 1.5 to 2. One can notice that going from 40 m 
to 60 m does not increase the signal-like events but in- 
creases the background-like events. For a 40 m length of 
the decay tunnel, the increase of the radius improves the 
number of signal- like events by 50%, and the backgroud 
increase by 70% to 100%. Notice that the v^jv^ flux ra- 
tio is rather insensitive to the decay tunnel length. The 
feeling that L T — 40 m and Rt = 2 m is a good signal 
over background compromise is confirmed by sensitivity 
quantitative studies reported in section 

Looking at the evolution of flux with respect to 
the beam energy, one notices that a maximum is reached 
around 4.5 GeV. This is due to the competition between 
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Table 5. Integral of the total flux of the different species with different settings. The f M an d fluxes are expressed in 
10 13 /100m 2 /y unit while the v e and z? e fluxes are expressed in 10 11 /100m 2 /y unit. The positive focusing and negative focusing 
are distinguished by a (+) sign and a (— ) sign, respectively. The settings used correspond to different values of L T and Rt, 
the length and radius of the decay tunnel. Setting (1) means Lt = 10 m and Rt = 1 m. Setting (2) is the default option and 
means Lt = 20 m and Rt = 1 m. Setting (3) means Lt = 20 m and Rt = 1.5 m. Setting (4) means Lt = 30 m and Rt = 1 m. 
Setting (5) means Lt = 40 m and Rt = 1 m. Setting (6) means Lt = 40 m and Rt = 1.5 m. Setting (7) means Lt = 40 m 
and Rt = 2 m. Setting (8) means Lt = 60 m and Rt = 1 m, and finally, setting (9) means Lt = 60 m and Rt = 1.5 m. 



Settings u e f M v e 









+ 




+ 




+ 




+ 




(1) 


2.2 


GeV 


5.5 


0.4 


1.7 


0.2 


0.3 


4.3 


0.1 


0.8 


(1) 


3.5 


GeV 


7.7 


0.7 


2.6 


0.6 


0.6 


6.6 


0.3 


1.3 


(1) 


4.5 


GeV 


7.1 


1.0 


2.8 


0.9 


0.5 


5.2 


0.3 


1.1 


(1) 


6.5 


GeV 


8.3 


1.2 


4.7 


1.9 


0.8 


5.6 


0.9 


1.8 


(1) 


8.0 


GeV 


7.7 


1.2 


5.1 


2.2 


0.9 


5.6 


1.1 


2.1 


(2) 


2.2 


GeV 


7.6 


0.4 


3.2 


0.2 


0.3 


5.8 


0.1 


1.6 


(2) 


3.5 


GeV 


10.0 


0.9 


4.4 


0.6 


0.7 


8.5 


0.3 


2.2 


(2) 


4.5 


GeV 


10.9 


1.1 


5.1 


1.0 


0.7 


6.7 


0.4 


1.8 


(2) 


6.5 


GeV 


10.4 


1.4 


6.4 


2.0 


1.0 


7.1 


0.9 


2.5 


(2) 


8.0 


GeV 


9.7 


1.5 


6.7 


2.3 


1.2 


7.1 


1.1 


2.8 


(3) 


2.2 


GeV 


9.0 


0.6 


4.4 


0.4 


0.4 


6.7 


0.2 


2.2 


(3) 


4.5 


GeV 


13.2 


1.5 


6.9 


1.4 


0.9 


8.1 


0.6 


2.7 


(4) 


3.5 


GeV 


10.9 


0.9 


5.7 


0.7 


0.7 


9.4 


0.3 


2.9 


(4) 


4.5 


GeV 


11.6 


1.2 


6.3 


1.0 


0.7 


7.1 


0.4 


2.3 


(5) 


2.2 


GeV 


8.9 


0.5 


5.1 


0.3 
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the cross section rise with respect to the energy and the 
decrease of the number of p.o.t due to the constant SPL 
power (4 MW). 

6 Sensitivity computation ingredients 

The sensitivity to #13 and Sep is computed for a — > v e 
appearance experiment. An analysis program described in 
reference ;37J has been used for such sensitivity computa- 
tion. See tableUOfor the default user parameter values used 
in this paper. We just remind here some key points of the 
program. 

It includes a full 3-flavors oscillation probability com- 
putation with matter effects, but no ambiguities are taken 
into account. This latest point may be revisited in a fu- 
ture work using reference J3$j. Concerning the background 
events, the v e jv e from the beam, the v^e~ elastic scatter- 
ing process, the ir° production as well as the /i/e misiden- 
tification are taken into account. The cross-sections from 
the NUANCE program are used PSj- The systematics er- 
ror on the total v e and v & fluxes determination is a user 
parameter and we have used the 2% value considered as a 



final goal, but also 5% and 10% |37j . The detector consid- 
ered for definitiveness is similar to the UNO detector, i.e. 
a 440 kt fiducial water Cerenkov detector |2()j . It is located 
at L = 130 km from CERN, in the foreseen new Frejus 
laboratory [21]. It is worth mentioning that if one wants 
to evaluate the influence of L on the sensitivity, it would 
mean a re-optimization of the horns for each L envisaged 
(see section Q}. The running time scenario has been fixed 
either by focusing positive particles during 5 years, either 
by focusing positive particles during 1 (or 2) year(s) fol- 
lowed by focusing negative particles during 4 (or 8) years. 

7 Results 

7.1 The positive only focusing scenario 

The 0i3 and S C p sensitivities are computed with 6> 13 = 
0° and Sep = 0° if not explicitly mentioned. It is worth 
stressing that the default parameters of table are used 
if not contrary mentioned, in particular, the decay tunnel 
geometry parameters (Lt — 20 m and Rt = 1 m), and 
the horn design to generate a 260 MeV neutrino beam. 
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Fig. 6. Comparison between the probability method, ( ) curve, and the full GEANT simulation method, ( ) curve, for 

the Vy, from ix + flux (a) and the P M from 7r _ flux (b). The horns are set to focus positive particles. It should be stressed that 
the full GEANT simulation has taken roughly 13 times more CPU time than the probability method with the same number of 
protons on target, and the later simulation is able to produce as well the v e and v e fluxes contrary to the former simulation. 



Table 6. Default parameters used to compute the sensitivity 
curves |3Z|- The quoted errors in parenthesis for the (12) and 
the (23) parameters (absolute value for the masse square differ- 
ences and relative value for the angles) are coming respectively 
from the up to date combined Solar and KamLAND results |4flj 
and from a 200 ktons-years SPL desappearance exposure |25j . 



Am\ 2 = 8.2(0.5) x 10" 


5 eV 2 sin 2 20i2 = 0.82(9%) 


Amis = 2.5(0.1) x 10" 


3 eV 2 sin 2 2623 = 1.0(1%) 


L T = 20 m 


Rt = 1 m 


M = 440 kT 


Csyst = 2% 


Horn/Reflector shapes 


to produce a 260 MeV neutrino beam 



Tabled presents the number of signal and background 
events for a 5 years positive focusing experiment, but with 
different beam energy settings. The significance parameter 
is defined in reference [HZj as 1 : 



yjN*? + (N*? xe syst ) 2 (1) 
with N*£ = N™ c + N% am + N othMkg 

and N°^ c the number of v e events due to oscillations, 
jybeam ^ ne num b er f background events coming from the 

1 Contrary to the definition of the significance of reference 
|37j . the systematical factor is applied to the total v e flux in 
agreement with the sensitivity contour computation. 



v e + D e contamination of the beam, N g the other 
kinds of background events and e sys t the systematical fac- 
tor. 

The contours at 90%, 95% and 99% CL of the 0i 3 sen- 
sitivity are presented in the (sin 2 2#i3, Z\m 2 3 ) plane on fig- 
ureEDfor 3.5 GeV proton beam kinetic energy. The com- 
parison between the contours at 90% CL with 2.2 GeV, 
3.5 GeV, 4.5 GeV and 8 GeV beam energies is shown on 
fignrelTll One notices in this scenario a better performance 
reached with a 4.5 GeV energy beam as a confirmation 
of significance parameter value. But, in fact there is not 
much visual difference between a sensitivity obtained with 
3.5 GeV and 4.5 GeV, even if one should keep in mind that 
kaon production models are different at these two energies 
(see section EJ. These two energy settings have been stud- 
ied with different decay tunnel geometry and results are 
reported on table EH One notices that similar results can 
be reached with a 3.5 GeV beam, compared to a 4.5 GeV 
beam. 

Quantitative studies of the minimum sin 2 26*i 3 with re- 
spect to the kinetic beam energy Ek(proton), and the de- 
cay length Lt, and the systematics e syst are presented in 
tables 03 and El One notices that for e syst = 5% there is 
no difference between a 3.5 GeV and a 4.5 GeV beam. 

We have also considered the 3.5 GeV and 4.5 GeV beam 
energies with the tunnel geometry parameters Lt — 40 m 
and Rt = 2 m, and the horn design producing a 350 MeV 
neutrino beam (see section QJ. In table are reported nu- 
merical values, and on figure El are shown the 90% CL 
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Fig. 7. Neutrino fluxes, 100 km from the target and with the horns focusing the positive particles. The fluxes are computed 
for a SPL proton beam of 2.2 GeV (4 MW), a decay tunnel with a length of 20 m and a radius of 1 m. The top left panel 
contains the f M fluxes, and the top right panel shows the fluxes. The bottom left panel presents the v a fluxes while the 

bottom right panel displays the i> e fluxes. The ( ) curve is the contribution from primary pions and the daughter muons, 

and from primary muons. The ( ) curve is the contribution from the charged kaon decay chain, and the ( ) curve is 

the contribution from the K° decay chain. An insert has been added to the plots to hight light when needed the contribution 
of charged and neutral kaons. 




1.2 1.4 
E (GeV) 



■ x10 7 



v„ flux 




1.2 1.4 
E (GeV) 



« 60 

0) 
>. 

S 50 

0) 

J 40 

I 30 
o 

? 20 
10 
0. 



X10 1 



v u flux 













x10 8 




50 










40 










30 










1 20 










1 10 










\ 0.2 0.4 0.6 0.8 1 1.21.4 



0.2 0.4 0.6 0.8 



V„ flux 



1.2 1.4 
E (GeV) 




1.2 1.4 
E (GeV) 



Fig. 8. Same legend as for figure but the horns are focusing negative particles. 
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Fig. 9. Same legend as for figure but for proton beam kinetic energy of 3.5 GeV (4 MW). 
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Fig. 10. Same legend as for figure |^1 but for proton beam kinetic energy of 3.5 GeV (4 MW). 
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Fig. 11. Same legend as for figure but for proton beam kinetic energy of 8 GeV (4 MW). 
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Fig. 12. Same legend as for figure |^1 but for proton beam kinetic energy of 8 GeV (4 MW). 
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Table 7. Number of events for 5 years positive focusing scenario with default parameters of table El Other backgrounds are 7T°, 
f M -elast., /i/e-missld. The significance parameter is defined by equation^ 





2.2 GeV 


3.5 GeV 


4.5 GeV 


6.5 GeV 


8 GeV 


non oscillated 


36917 


60969 


73202 


78024 


76068 


oscillated u e 


43 


60 


64 


61 


56 


beam u e 


165 


222 


242 


288 


299 


other background 


70 


105 


127 


148 


152 


Significance 


1.88 


2.16 


2.17 


1.87 


1.69 



Table 8. Minimum sin 2 2^13 x 10 3 in the (sin 2 2#i3, Z\m| 3 ) plane observable at 90% CL computed for different decay tunnel 
length (Lt) and kinetic beam energy (E k (proton)) and 5 year of positive focusing. Other parameters are fixed to default values 
(table 0. 
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Fig. 13. Sensitivity contours obtained with a SPL energy of 
3.5 GeV and default parameters of table El In particular, it is 
reminded that the tunnel geometry parameters are Lt = 20 m 

and Rt = 1 m. ( ), ( ) and ( ) curves stand for 

90%, 95% and 99% confidence level, respectively. 




sin 2 2i? 13 

Fig. 14. Comparison of 90% CL sensitivity contours obtained 

with SPL energies of 2.2 GeV ( ), 3.5 GeV ( ), 

4.5 GeV ( ) and 8 GeV ( ) and default parameters of 

table El In particular, it is reminded that the tunnel geometry 
parameters are Lt = 20 m and Rt = 1 m. 



sensitivity contours. With the 350 MeV neutrino beam, 
one can expect a 16% improvment with respect to the 
260 MeV neutrino beam for the same decay geometry. 
One also notices that there is marginal gain to increase 
the beam energy from 3.5 GeV to 4.5 GeV, as already 
mentioned. 

As well, there are variations on the minimum sin 2 2#i3 
value that may be reached in a — ► v e experiment which 
are due to the sign(Z\m2 3 ) ambiguity and the Sep value. 



On table El are presented these kinds of variations. Other 
ambiguities coming from the sign(tan(26>23)) ignorance also 
occur as studied in reference From figure 9 of this 
reference, we estimate a 30% effect on sin 2 (26*13) sensivity 
due to these ambiguities. 
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Table 9. Minimum sin 2 2^13 x 10 3 in the (sin 2 2#i3, Z\m 2 3 ) plane observable at 90% CL computed for different decay tunnel 
length (Lt) and radius (Rt) for the 3.5 GeV and 4.5 GeV scenarios and 5 years of positive focusing. Other parameters are fixed 
to default values (table EJ- Settings in parenthesis are identical to those of table except that the setting (7b) corresponds to 
the tunnel geometry of setting (7) but the horn geometry producing a 350 MeV neutrino beam is used. We remind that the 
setting (2) is the default one, and the settings (7) and (7b) correspond to Lt = 40 m and Rt — 2 m. 



setting 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (7b) 


(8) 


(9) 


3.5 GeV 
4.5 GeV 


0.92 
1.04 


0.92 
0.89 


0.83 
0.82 


0.98 
0.94 


1.00 
0.99 


0.93 
0.92 


0.91 0.76 
0.87 0.71 


1.05 
1.03 


1.01 
1.00 




sin 2 2^ 13 

Fig. 15. Comparison of 90% CL sensitivity contours obtained 
with SPL energies of 3.5 GeV or 4.5 GeV, and either a 260 MeV 
(default) neutrino beam or a 350 MeV neutrino beam. The 
tunnel geometry parameters are Lt = 40 m and Rt = 2 m. 
The ( — ■ — ) curve corresponds to a 350 MeV/4.5 GeV (neu- 
trino beam/SPL beam energy) setting; the ( ) curve cor- 
responds to a 350 MeV/3.5 GeV setting; the ( ) curve cor- 
responds to a 260 MeV/4.5 GeV setting and the ( ) curve 

corresponds to a 260 MeV/3.5 GeV setting. 

7.2 Mixed positive/negative focusing scenario 

The combined sin 2 29i3 and Sep sensitivity for the 5 years 
positive focusing scenario and the default parameters of 
table El is presented on figure ITfiT a). The results obtained 
with a 3.5 GeV and 4.5 GeV SPL beam are similar and 
better than with the other energy settings. On figureEJa) 
the results obtained with a 260 MeV neutrino beam and 
a 350 MeV neutrino beam are presented with a 40 m 
long, 2 m radius decay tunnel. With the 350 MeV neu- 
trino beam, one can reached better sensitivity results in 
the range \Scp\ < 120°, and comparatively the gain ob- 
tained when switching from a 3.5 GeV proton beam to a 
4.5 GeV proton beam is marginal. To improve the Scp- 
independent limit on sin 2 26*i3, especially around Sep = 



Table 10. Minimum sin 2 29 1S x 10 3 in the (sin 2 26 13 , Am\ 3 ) 
plane observable at 90% CL computed for different level of sys- 
tematics (e a yat) and kinetic beam energy (Ek(proton)) and 5 
years of positive focusing. Other parameters are fixed to default 
values (table 0. 





2.2 GeV 


3.5 GeV 


4.5 GeV 


6.5 GeV 


8 GeV 


2% 


1.16 


0.92 


0.89 


1.01 


1.12 


5% 


1.48 


1.25 


1.25 


1.48 


1.64 


10% 


2.40 


2.14 


2.21 


2.72 


3.09 



Table 11. Minimum sin 2 26>i 3 x 10 3 in the (sin 2 20i 3 , Amis) 
plane observable at 90% CL computed for a 2.2 GeV kinetic 
energy proton beam, and for different values of sign(Z\m| 3 ) 
and Sep and 5 years of positive focusing. Other parameters 
are fixed to default values (table EJ- 



-180° 


-90° 


0° 


90° 


180° 


+ 1.40 


0.43 


1.16 


11.48 


1.40 


- 1.45 


11.75 


1.11 


0.43 


1.45 



90°, one may envisage a combination of 2 years with pos- 
itive focusing and 8 years negative focusing as in refer- 
ences [T7ll25ll26| . The comparison of the results obtained 
with different SPL beam energies on the combined sensi- 
tivity contours are presented in figure lTfiT h). Quantitative 
results with this kind of mixed focusing scenario are re- 
ported table El O ne generally gets 10% to 20% better 
limit on sin 2 2#i 3 independently of Sep with a 3.5 GeV 
kinetic energy beam compared to a 2.2 GeV beam. Dou- 
bling the length and the radius of the decay tunnel allows 
to reach a 10% better limit. 

On figure ITTT b) are presented the results considering 
the effects of a 350 MeV neutrino beam obtained either 
with a 3.5 GeV proton beam or a 4.5 GeV proton beam 
compared to a 260 MeV neutrino beam obtained with a 
4.5 GeV proton beam. The tunnel geometry parameters 
are Lt = 40 m and Rt = 2 m (other tunnel geometry 
have been studied but the results are worse and so are 
not reported). Except in the region \Scp\ > 150°, the re- 
sults obtained with the 350 MeV neutrino beam (3.5 GeV 
proton beam) are somewhat better, even if a 11% im- 
provement of the c> C p-independent sin 2 2# 13 limit can be 
reached with the 260 MeV neutrino beam obtained with 
the 3.5 GeV proton beam. 
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Fig. 16. 90% sensitivity contours obtained with SPL beam 

energy of 2.2 GeV ( ), 3.5 GeV ( ), 4.5 GeV ( ) 

and 8 GeV ( ) at 90% CL. Default parameters of tabled 

are used either with a 5 years positive focusing scenario (a) 
or a mixed scenario of 2 years positive focusing and 8 years 
of negative focusing (b). In particular, it is reminded that the 
tunnel geometry parameters are Lt = 20 m and Rt = 1 m. 

8 Summary and outlook 
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Fig. 17. 90% CL sensitivity contours obtained with the de- 
cay tunnel geometry parameters Lt = 40 m and Rt = 2 m 
and different SPL beam energies (3.5 GeV or 4.5 GeV) and 
different horn designs (260 MeV or 350 MeV neutrino beams): 

( ) curve for a 350 MeV/4.5 GeV setting, ( ) curve for a 

350 MeV/3.5 GeV setting, ( ) curve for a 260 MeV/3.5 GeV 

setting. Other default parameters of table are used either 
with a 5 years positive focusing scenario (a) or a mixed scenario 
of 2 years positive focusing and 8 years of negative focusing (b). 



A complete chain of simulation has been set up for the 
SPL-Frejus project. The neutrino production has been ex- 
tended to the kaon decay contribution, which is important 
to test SPL energy scenario above 2.2 GeV. 

The beam line optimization has been performed us- 
ing the sensitivity to sin 2 2613 and Sep- The shape of the 
focusing system has been updated to obtain a neutrino 
beam energy around 260 MeV or 350 MeV. 

In a positive only focusing scenario, the best limit on 
sin 2 26>i3 is 0.71 x 10~ 3 (90% CL, S C p = 0), with a 4.5 GeV 
beam energy and a 40 m long, 2 m radius decay tunnel, 
and a beam energy around 350 MeV. However, the 3.5 GeV 
beam may also obtain rather similar limit with 0.76 x 
10~ 3 (90% CL) with the same tunnel parameters. The 
Sep independant sin 2 26^3 sensitivity is limited to « 10 -2 
due to the 8c p ~ 90° region. 

In a mixed focusing scenario, the best limit on sin 2 26 13 
independent of S C p is 2.02 x 10~ 3 (90% CL) obtained 
with a 3.5 GeV beam energy, a 40 m long, 2 m radius 
decay tunnel and a beam energy around 260 MeV. But 
for \Scp\ < 150°, the 350 MeV neutrino beam is better, 
keeping the primary proton energy at 3.5 GeV and with 
the same decay tunnel parameter. 

The comparison of the optimization presented in this 
paper with the results obtained by other projects is dis- 
played on figure El It presents the 5 years positive focus- 



ing scenario, and two versions of a mixed scenario using 
positive and negative focusing: one scenario duration is 5 
years in total and the other one is 10 years running in 
total and has been used in the previous section. It shows 
the complementarity of the SPL-Frejus project with the 
beta beam-Frejus project. Especially when considering the 
sensitivity to sin 2 2813 for Sep < 0. 

The authors think that the present study may be ex- 
tended in many respects. The beam line simulation part 
may be performed with a single simulator as FLUKA (or 
GEANT4 [45. a for comparison). Other targets may be en- 
visaged (tantalum, carbon) as well as other detector types 
as a Large Liquid Argon detector The baseline length 
may also be revisited as well as the off axis option. The 
sensitivity analysis may be deeper investigated using the 
complete set of possible ambiguities as in reference [T7] . 
and the 6*13 or Sep measurement accuracy with new beam 
energy scenario may be investigated too. 



The authors would like to thank M. Mezzetto for expressing 
his interest since the early stage of this work and for providing 
us with his sensitivity computation program. Also the authors 
thank S. Gilardoni for fruitful discussions. 
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Table 12. Minimum sin 2 2#i3 x 10 3 observable at 90% CL computed for the worse Scp case, and for different decay tunnel 
length (Lt) and radius (Rt) for the 3.5 GeV and 4.5 GeV scenarios and 2 years of positive focusing plus 8 years of negative 
focusing. Other parameters are fixed to default values (table EJ- Settings are identical to those of table 



setting 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(7b) 


(8) 


(9) 


2.2 GeV 




2.52 






2.58 




2.30 








3.5 GeV 


2.34 


2.22 


2.10 


2.13 


2.09 


2.08 


2.02 


2.28 


2.16 


2.09 


4.5 GeV 


2.91 


2.60 


2.43 


2.48 


2.52 


2.39 


2.34 


2.55 


2.53 


2.47 
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Fig. 18. 90%CL sensitivity contours labeled by the project or 
experiment involved. The "CHOOZ excluded" dashed curve 
comes from the exclusion obtained from reference flj] with 
Am 2 = Zim 2 tm ; in the same conditions is given the sen- 
sitivity foreseen for the "Double-CHOOZ" project [HI- The 
"CNGS combined" has been obtained combining the results 
form OPERA and ICARUS gT]. The T2K contour has been 
derived from reference [42 j . The BNL contour has been ob- 
tained from reference |43j . The "Beta Beam" contour has been 
computed with 5 years running with both u e and u e neutrino 
beams in an appearance mode, while the dashed "Beta Beam 
disappearance" has been obtained as if the f3 beam were anal- 
ysed like a reactor experiment with 1% systematic error |44j . 
The "SPL 5y" and "SPL 2y+8y" and "SPL ly+4y" curves 
have been obtained from the optimisation described in this 
paper ("5y": positive only focusing scenario; "ly+4y": 1 year 
of positive focusing and 4 years of negative focusing scenario; 
"2y+8y": 2 years of positive focusing and 8 years of negative 
focusing scenario) using a 3.5 GeV beam and a decay tunnel 
of 40 m length, and 2 m radius. 



A Decay probability computations 

This appendix contains the probability formulas and the 
algorithms used in the flux computation (see section I57T)) . 



A.l Pion neutrino probability computation 

Pions decay only as ir + — > /i + + or 7r~ — > fx" + and 
the neutrinos are emitted isotropically in the pion rest 
frame, with an energy of about 30 MeV given by the 2- 
body decay kinematics. Applying a Lorentz boost knowing 




Fig. 19. Pion decay in the tunnel frame. To reach the detector, 
5 = — a is needed. 



the pion momentum and direction, it is possible to com- 
pute the probability to reach for the neutrinos the detec- 
tor. Only neutrinos parallel to the beam axis are supposed 
to pass through the detector fiducial area, and therefore, 
the neutrinos must be emitted by the pion with an angle 
opposite to the angle between the pion and the beam axis 
(see figure HSjl . This gives: 



P* = 



1 A 



l-(3 2 



47T L 2 (Pcosa - 1) 



(2) 



where (3 is the velocity of the pion in the tunnel frame, A 
is the fiducial detector surface, L the distance between the 
neutrino source and the detector, and a the angle between 
the pion direction and the beam axis in the laboratory 
frame. 



A. 2 Muon neutrino probability computation 

Muons decay only as fi + — > e + + v e + or (x~~ — > e~ + v e + 
v^, and will produce background events. The mean decay 
length of the muons is 2 km, therefore, most of them do 
not decay in the tunnel. This induces a lack of statistics 
to estimate the corresponding level of background. This 
problem has been solved using each muon appearing in 
the simulation in the following steps: 

1. the probability for the muon to decay into the tunnel 
has been computed using a straight line propagation; 

2. the available energy for the neutrino in the tunnel 
frame has been divided in 20 MeV energy bins; 

3. one v e and one have been simulated in each of the 
energy bins (step 2). Then, the probability to reach 
the detector has been computed, and multiplied by 
the probability computed at step (1). 

After the probability computation, the non useful muon 
is discarded by GEANT to gain in CPU time. 
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Table 13. Flux function in the muon rest frame |47j . 



fo(x) 



2x 2 (3~2x) 2x 2 (l~2x) 



12x^(1 -x) I2x 2 (l-x) 



Table 14. Charged and neutral kaon decay channels |49j . 



A" 





63.51% 


7v~e + v e 


19.35% 7T + 7T- 


68.61% 


n ± n° 


21.17% 


■K + e~v e 


19.35% ttV 


31.39% 




5.59% 




13.5% 






4.82% 


7T + >~P M 


13.5% 






3.18% 


7T 7T 7T 


21.5% 




7Y 7T 7Y 


1.73% 




12.38% 





The probability for the muon neutrino and the electron 
neutrino to be emitted parallel to the beam axis is 36J: 



dE v 



1 A 2 



1 



4tt L 2 to m7ai (1 + /^cosfl*) 



order to obtain statistically significant results. The mul- 
tiplicity of decay chann els m akes impossible the method 
used for the muon case l|A.2j) . The method chosen for the 
good compromise between the gain in CPU and the statis- 
tical uncertainty of the results, is to duplicate many times 
each kaon exiting the target. The number of duplication 
varies between 10 and 300. It depends on the initial kaon 
rate and therefore on the beam energy. 

All the kaons daughter particles are tracked by GEANT 
until they decay. Three different types of daughter parti- 
cles are identified in the kaon decays. The first type cor- 
responds to primary neutrinos, the second type concerns 
charged pions and muons, and the neutral pions are left 
for the last type. 

In the — > ^v^ip^ decay modes, the computation 
of the probability for a neutrino to reach the detector is 
the same than the 2-body decay formula used to in the 
pion decay (equation^, where (3 is now the kaon velocity, 
and a the angle of the kaon with respect to the beam axis. 

When a neutrino is produced by a kaon 3-body decay, 
the probability to reach the detector is computed using a 
pure phase space formula. It yields: 



Hi 



dV K 



1 A 



1 



(/3 M cosp- 1) : 



[/ o (s)TJ#/i(aOcos0*] (3) dE v in L 2 m K 



uin - mi 
1 



where /? M and 7 M are the velocity and the Lorentz boost of 
the muon in the tunnel frame, 9* is the angle with respect 
to the beam axis of the muon in the muon rest frame, 
p is the corresponding angle in the tunnel frame. As the 
pion case, this angle appears because the neutrino must 
be parallel to the beam axis. 77^ is the muon longitudinal 
polarization, the parameter x is defined as x — 2E*/m^ 
where E* is the neutrino energy in the muon rest frame, 
and the functions fo(x) and fi(x) coming from the matrix 
element of the muon decays are given in tabled The sign 
in front of 77^ in equation's (— ) for the fi + decays and 
(+) for the fx~ decays, respectively. 

Muon polarization is computed using the conservation 
of the transverse component of the velocity four-vector 
7(1, (3) between the muon rest frame (where the polar- 
ization is computed) and the pion rest frame, where the 
muon helicity is —1, due to the parity non conservation. 
It yields gHj: 



0: 



7x(l+/3ifCos<9*) (/3 K cosS-1) 2 



(5) 



7aA 



sin6>* and II . 



(4) 



where 7^, /?„•, 7^, and (3^ are the Lorentz boost and veloc- 
ity of the pion and of the muon in the tunnel frame, and 
9* the angle with respect to the beam axis of the muon in 
the pion rest frame. 



A. 3 The treatment of the kaons 

Contrary to pions and muons, kaons have many decay 
channels. They are summarized in table ITU 

There is a very small amount of kaons produced (sec- 
tion EJ, and this number has been artificially increased in 



where mx is the kaon mass (charged or neutral), is 
the pion mass (it mass in K ± decays and n ± mass in 
Kl decays), and to/ is the mass of the lepton associated 
with the neutrino. The (3k and jk are the velocity and 
the Lorentz boost of the kaon, 9* is the angle between the 
neutrino direction and the kaon direction, in the kaon rest 
frame. Finally, S is the angle between the kaon direction 
and the beam axis in the tunnel frame. 

When a tt^ 1 is produced in the kaon decay chain, it is 
tracked by GEANT until it decays, and the probability 
of equation El is applied to the produced neutrino. In case 
of a muon, it is treated as explained in IA.2I The muon 
polarization is computed this time using the kaon decay 
informations. Finally, when a ir° is produced, as it cannot 
create neutrinos, it is simply discarded. 
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